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Abstract 

A systematic study of the pairing-correlations derived from various particle- 
number projection methods is performed in an exactly soluble cranked- 
deformed shell model Hamiltonian. It is shown that most of the approximate 
particle-number projection methods including the method of Lipkin-Nogami, 
which is used quite extensively in nuclear structure studies, breakdown in 
the weak pairing limit. The results obtained from the recently formulated 
number-projected Hartree-Fock-Bogoliubov (PHFB) equations, on the other 
hand, are in complete agreement with the exact solutions of the model Hamil- 
tonian. The pairing-energy calculated from PHFB method is shown to be 
finite in all the studied limits. More importantly, the numerical work in- 
volved in the solution of the PHFB equations appears to be comparable to 
the solution of the bare HFB equations. 

PACS numbers : 21.60.Cs, 21.10.Hw, 21.10.Ky, 27.50. +e 
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I. INTRODUCTION 



Pairing-correlations play a central role in describing the properties of atomic nuclei ■ 
The bulk of the nuclear phenomena, for instance, the odd-even mass differences and moments 
of inertia of deformed nuclei can be described by considering nucleons to be in a superfluid 
paired phase. The recent discovery of halo nuclei have also elucidated the importance of 
pairing-correlations in the formation of such structures In the traditional nuclear 

structure models, the pairing-correlations are treated in Bardeen-Cooper-Schrieffer (BCS) or 
Hartree-Fock-Bogoliubov (HFB) approximations 0. These approaches use a simple ansatz 
of a product-state of quasiparticles for the ground-state wavefunction. The product ansatz 
leads to a set of non-linear equations in the BCS theory or to the diagonalization of a matrix 
in the case of HFB approach. There exist now very fast numerical methods to find accurate 
solutions of these problems. The major advantage in the mean-field models is that the 
number of non-linear equations or the dimensionality of the matrix to be diagonalised is 
quite small as compared to the exact analysis, which is impossible to apply in a realistic 
configuration space of medium and heavy mass nuclei. 

Nevertheless, the product ansatz of the mean-field approach breaks the symmetries which 
the orginal many-body Hamiltonian obeys. For instance, in the case of BCS theory, the 
product-state violates the gauge-symmetery associated with the particle-number. The BCS 
product-state of quasiparticles does not have a well defined particle-number. The most 
undesirable feature associated with the broken symmetry is that it leads to sharp phase 
transitions. These phase transitions may be reasonable for a system with a very large 
number of particles, for example, a metallic superconductor. But for a finite mesoscopic 
system these phase transitions are washed out by fluctuations and are not observed in the 
experimental data or in the exact study of toy models. The sharp transition is an artifact 
of the simple ansatz of the product-state, which does not contain fluctuations. 

The symmetry restoration in the mean-field studies is a major challenge to the nuclear 
structure models. Although, the projection methods which restore the broken symmetries 
have been known for more than thirty years, but have been studied only in simple model cases 
due to the numerical difficulties. In most of the realistic studies, approximate projection 
methods have been employed 0. Only recently, the gradient method in the projected- 
energy surface has been employed in realistic calculations with Gogny forces There has 
been an unsolved problem whether a simple HFB like equations can be obtained with the 
projected-energy functional. This would allow the usage of fast numerical algorithms for 
diagonalisation of matrices as is done in the ordinary HFB theory. This problem has been 
recently addressed and it has been shown that the variation of an arbitrary energy-functional, 
which is completely expressed in terms of the HFB density matrix p and the pairing tensor 
K, results in the HFB like equation It is noted that the projected-energy functional can 
also be expressed in terms of the HFB density matrices p and k and, therefore, one obtains 
HFB like equations with modified expressions for the pairing-potential and the Hartree-Fock 
field. 

In the present work, we shall limit ourselves to the restoration of the gauge symmetry 
associated with the particle-number projection, which is important in the discussion of the 
pairing-correlations. Several approximate particle-number projection methods have been 
developed and the most popular of these is the method of Lipkin and Nogami 0. The 
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purpose of the present work is to perform a systematic study of these various approximate 
projection methods and to compare them with the recently developed number-projected 
HFB (PHFB) approach. It is shown that most of the approximate methods including the 
Lipkin-Nogami method breakdown in the weak pairing limit. The PHFB method, on the 
other hand, gives finite pair-correlations in all the cases studied. The paper is organised as 
follows : a brief survey of the the particle-number projection methods is given in the next 
section and in the appendix. The model Hamitonian used in the present work is expressed 
in section III, the numerical results are discussed in section IV and finally the present work 
is summarised in section V. 



II. PROJECTION METHODS 

In this section, we shall briefly describe the approximate particle-number projection 
methods used in the literature and the exact number-projected HFB method proposed re- 
cently P). Although, the approximate methods are quite standard and the details can be 
found, for example, in the text book [Q. But for completeness, the essential elements of 
these methods are discussed in the following subsections. 



A. HFB method 



The many-body Hamiltonian of a system of fermions is usually expressed in terms of a 

I, ■ ■ ■ , Cj^; Ci, . . . , %^): 



set of annihilation and creation operators (c, c^) = (c^, . . . , Cj^; 



H = N p c + — \ 7; c c f 1 ^ 

/ J ^nin2 ni^n2 ' a / ^ ^nin2n3ri4 ni n2 n4^n3' \ / 

where the anti-symmetrized two-body matrix-element is defined by 

^nin2n3n4 = ("'l"'2 1 1^ |n3n4 - Uim) . (2) 

The operators annihilate the bare fermion vacuum 

c„|-) = 0, (3) 

for all n and provide a basis in the present analysis. 

The most basic of the mean-field approaches to describe the pairing-correlations is the 
HFB method. In this method, the ground-state wavefunction is a product-state of quasi- 
particle operators {a,a^) = {ai, . . . ,a^;ai, . . . ,a\^). These quasiparticle operators are 
connected to the original particle operators by a linear transformation 

«fc = II {KkCn + K*fc4) , (4) 

n 

«I = E (K^feC„ + f/„fe4) . (5) 



This transformation can be rewritten in the matrix form as 
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The quasiparticle operators need to satisfy the same fermion commutation relations as the 
original operators in order to preserve the Fermi-Dirac statistics and, therefore, the trans- 
formation matrix is required to be unitary 

W^W = WW^ = /, (7) 

which leads to following relations among the coefficients Unk and Vnk, 

U^U + = I, UU^ + V*V^ = I, (8) 
U'^V + V^U = 0, UV^ + V*U'^ = 0. (9) 

The quasiparticle operators annihilate the quasiparticle vacuum |$), defined by 

«fc|*) = 0, (10) 

for all k. In mean-field theory, it represents an approximation to the ground-state of a 
system and turns out to be a generalized Slater-determinant 0. Since the quasiparticle 
transformation mixes creation and annihilation operators, |$) does not correspond to a 
wavefunction with good particle-number. In HFB theory, we have two kinds of densities, 
the normal density p, and the pairing-tensor k, defined as 

Pnn' = ($|4c„|<l>), Knn, = ($ | C,,C„ | $) . (11) 

These can be expressed in terms of the HFB coefficients as 

p = V*V^, K = V*U^ = -UVl (12) 
The HFB energy is given by 

where the HFB fields are defined as 

(14) 

n2n4. 



The variation of the HFB energy-functional (|T3D results in the standard HFB equations : 
where 



^(1]=eAI]. (16) 



^= ('^-A* ^ -(e + r)* + A)- 

As already mentioned, the HFB wavefunction does not have a well defined particle-number 
and the Lagrangian parameter (A) is introduced to have the correct particle-number on 
the average. The diagonalisation of the HFB matrix ( p!?] ) gives rise to the quasiparticle 
energies {Ei) and the excitation spectrum of the many-body system is constructed from 
these quasiparticle energies. 
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B. Variation after projection - Number-projected (PHFB) equations 



It has been shown recently that the variation after projection leads to the HFB like 
equations with modified expressions for the pairing-field and the Hartree-Fock potential. We 
have also obtained these equations using the method of statistical averages and the derivation 
is given in the appendix. In this section, we shall present the relevant expressions of the 
projected-fields which are used in the numerical study. The project ed-energy functional is 
given by 



where p and n are the unprojected densities given in Eq. |T2| and the particle-number 
projection operator is defined as 

pN = l_ /rf0e*0{^-^). (19) 

271 J 

It has been demonstrated in ref. |P that the variation of the projected-energy results in the 
HFB equations of the same structure as in ( p!7D except that the expressions for the fields are 
different. The projected HFB equation is given by 



where 



^ ~ y — (A^)* — (e^)* — (r^)* — (A^)* j ' ^^''"^ 
The number-projected expressions for the fields are now given by 



^JdcP y{<P) {F(0)Tr[ep(0)] + [1 - 2^e-^<^ sin 0p(0)]eC 



+ h.c. (22) 

■ j d<P yi<P) (r(</))iTr[r(0)p(0)] + ^[1 - 2ze-*^sin#(0)]r(0)C(0)) 

+ h.c. (23) 

I d<P y{<P) (F(0)^Tr[A(0)7^*(0)] - 2^6"^^ sin C(0) A(0)r) 
+ h.c. (24) 



A^ = i / y{ct^)e-'"fC (</.) A(</.) - (..)^, (25) 



with 
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(0) = (0), (26) 



K^n2i.'P) = \Y. Vn2n3n4'^n3n4(<^)' (^7) 
n3n4 

K^nS^t^) = O <in2(0)^nin2n3n4) (28) 

nin2 

p(0) = C(0)p, (29) 

k(0) = C(0)k = kC^(0), (30) 

^(0) = e^'"^KC*{(P) = e^"^C\(P)K, (31) 

C(0) = e2^<^(l + p(e2**-l))"\ (32) 

X = ^ , (33) 



yi<^)-J^y ld9yi<P)-l, (34) 



F(0) = ie-*-^ sin C(0) - z / rf0'i/(0')e"*^' sin 0' C(0'). (35) 



and 



The term designated by in Eq. (p4D does not appear in the normal HFB formalism 
and it can be immediately shown that it vanishes for the gauge-angle, = 0. This term 
originates from the variation of the pairing-energy with respect to the density-matrix. In 
normal HFB theory, the pairing-energy depends only on the pairing-tensor, but the projected 
HFB pairing-energy also depends on the density-matrix through the norm-overlap. As a 
matter of fact, in a general case, the norm-overlap depends on both the density-matrix and 
the pairing-tensor (see appendix). But for the special case of number-projection, the term 
in the overlap-matrix, depending on the pairing-tensor, can be rewritten in terms of the 
density-matrix by using the HFB relation {p — = kk,^). Due to this transformation, the 
expression for in (25) has a very simple appearance and reduces to the familiar form in 
the canonical representation |l8|,[T0[|. 

It should be noted that as compared to the HFB case, the projected quasiparticle energies 
obtained from the diagonalisation of the projected HFB matrix Eq. (^) have no meaning. 
The only meaningful quantity is the total projected-energy, which is given by 

Etot = J rf0y(0)ff(0), (36) 

■ dM^) {Hspi^) + Hpf^i^) + Hp^m , (37) 



where. 



if,p(0) = Tr(ep(0)), (38) 
i7,,(0) = ^Tr(r(0)p(0)), (39) 

^Pp(0) = -^Tr(A(0)r(0)). (40) 
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In the present work, we shall discuss the pairing-energy (Epair) and the aligned angular- 
momentum (< Jj. >) quite extensively, since these give a measure of the pairing-correlations. 
The expressions for these quantities are given by 



Epair = J d(j)Hpp{(j)), (41) 

< J,. > = /rf0Tr(J,.p(0)). (42) 



As is clear from Eq. (|21|), the projected equations have exactly the same structure as 
that of normal HFB. Therefore, one can employ the existing HFB computer codes and 
only the expressions for the projected- fields need to be evaluated. The projected-fields now 
involve the integration over the gauge-angle, 0. This integration has been performed using 



the Gauss- Chebyshev quadrature method |11]. In this method, the integration over the 



gauge-angle is replaced by a summation. It can be shown |TT| that the optimal number of 
mesh-points in the summation which eliminates all the components having undesired particle 
numbers is given by 

M = max Qat, ^] - + 1, (43) 
where N is the number of particles and Q is the degeneracy of the model space. 



C. Projection after variation (PAV) 

In this projection method, the HFB mean-field solution is used to calculate the projected 
quantities. The normal variational calculations are carried out in the HFB formalism and 
the solution obtained is then used to calculate the projected quantities as given by Eqs. 



(^,^,|^). Therefore, this method is closely based on the HFB method and in the event 
that HFB solution depicts a transitional behaviour, it would also be reflected in the PAV 
results. 



D. Lipkin-Nogami (LN) method 

The Lipkin-Nogami (LN) method has been used quite extensively as an approximation 
to the particle-number projection approach. In this method, the HFB energy-function is 
modified to include the quadratic term in the Kamlah expansion of the projected-energy [Q. 
The HFB energy is replaced by 

Ehfb Ehfb - X2{AN^), (44) 
where AiV^ = (iV — {N))^ and the Lagrangian parameter A2 is given by |0 

'~ (AiV4) - {ANY - (AiV3)V(AiV2) ■ ^ ^ 

In the above equations, the expectation value of an operator is calculated with respect to 
the HFB vacuum, i.e., (O) = (<I>|0|<I>) and can be calculated by differentiating the norm 
and the Hamiltonian overlaps as a function of the gauge-angle, i.e.. 
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(46) 



1 dPH{(j)) 



\<t>=Oi 



(47) 



where x{(j)) and H[(j)) are defined in Eqs. ( P3| ) and (PBD, respectively. The expressions for 
the expectation values required in Eq. (^) are given, for example in ref. 0. We have 
rederived these expectation values using the expressions (|33D and (|36|) , which essentially 
involve differentiation of determinants and the matrix functions. The expectation values are 
given by 



{AN'') 
(AiV^) 
{AN^) 
(HAN) 

{AN'))) 



2Tr(x), 
4Tr(7x), 

3(AiV2)2 + 8Tr{x(l-6x)}, 
2Tr(/ix) -7^Tr(Afi:*7), 

4TT[{h-f + VHF{x})x], 

7^Tr[Afi;*(l - 8x)] - Tr[Vp{7K}(7K)*] 



where, 



X = p(i -p), 

7 = 1 - 2p, 
h = e + r, 

VHF{x}ij = ^VikjlXlj, 

kl 

Vp{'^K}ij = ]:J2'^ijkl{l'^)kl- 
^ kl 



(48) 
(49) 
(50) 
(51) 

(52) 



(53) 
(54) 
(55) 
(56) 

(57) 



TZ in Eq. ([52|) denotes the real part of the expression. In the LN method, the calculations 
proceed exactly as in HFB case except for the extra A2 term. At each HFB iteration, the 
A2 is calculated using the expression (^Sj) and this is then used in the solution of the HFB 
equations in the next iteration. 



III. MODEL HAMILTONIAN 

In the present work, we have carried out extensive numerical calculations using the 
model Hamiltonian consisting of a cranked deformed one-body term, h' and a scalar two- 
body delta- interaction, V2 [|12|. The deformed shell model Hamiltonian employed is given 

by 

H = + V,, (58) 

= hdef -UjJx + V2, (59) 

where. 
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V m 



and 



Vo 



LM 



(60) 



(61) 



with A\ 



'-LM 

matrix-element 



CjCj)LM 



and A 



LM 



(4 



LM) 



For the antisymmetric-normahzed two-body 



Ej ), we use the delta-interaction which for a single j-shell is given by |T 



Er 



-G 



(2j 



2{2L + 1) 



3 J L 

2 2 " 



(62) 



where the symbol [ ] denotes the Clebsch- Gordon coefficient. We use G = g J R^f'^dr as 
our energy unit and the deformation energy n is related to the deformation parameter f3. 
For the case of /iii/2 shell, k=2.4 approximately corresponds to (3 = 0.23 |T^. 

Although, the present model Hamiltonian is not very realistic and the results obtained 
cannot be compared with the experimental data. Nevertheless, it contains most of the 
essential components of a realistic Hamiltonian. The advantage in this model is that it can 
be solved exactly and it is possible to test the quality of an approximate scheme. We consider 
that it is quite appropriate to test the number-projection methods in a cranking model, since 
with increasing rotational frequency the pairing- correlations drop and the results become 
quite sensitive to the projection method. 



IV. RESULTS AND DISCUSSIONS 



The detailed numerical calculations have been carried out for six-particles in j = 11/2 
shell in order to investigate the validity of the particle-number projection methods pre- 
sented in section II. The calculations have been done with three different pairing-interactions, 
monopole [L=0 in Eq. (p^)], monopole plus quadupole [L=0 and 2 in Eq. (p2|)] and with 



full delta-interaction to check the sensitivity of the projection method on the underlying 
pairing-interaction. 

First, we would like to mention that the projected-HFB calculations almost completely 
agree with the exact solutions for all the cases studied and it is not possible to distinguish 
them in the figures. In all the figures, we present only the PHFB results and it should 
be kept in mind that these results coincide with the exact solution. In Fig. 1, the results 
of the total-energy {Etot) and the pairing-energy (Epair) are presented as a function of the 
strength of the pairing-interaction. For this set of calculations, we have employed the full 
delta-interaction. The HFB total-energy deviates from PHFB with increasing strength of 
the pairing-interaction. The results of total-energy with projection after variation appear to 
be in good agreement with the results of PHFB. The calculated total-energy in LN method 
interpolate between the PHFB and the HFB results. In general, it is noted from Fig. 1 that 
the total-energy is fairly reproduced by all the approximate projection methods. 

The critical test of a projection method is how well it reproduces the pairing-energy. In 
the lower panel of Fig. 1, the pairing-energy obtained from different projection methods is 
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presented. The major problem associated with the HFB theory is that the pairing-cncrgy 
shows a phase transition from finite pairing to zero pairing as is clearly evident from Fig. 1. 
The HFB pairing-energy is finite for pairing-strength, G greater than 0.4, but collapses for 
G less than or equal to 0.4. The PAV projected-energy has a similar behaviour as that of 
HFB energy and is expected since it uses the HFB mean-field solution. It should be noted 
that the PAV pairing-energy for certain values of the strength has a higher value than the 
HFB pairing-energy. 

The LN pairing-energy is somewhat between HFB and PHFB pairing- energy and has 
a non-zero value for G — 0.4. But for the lower values of the pairing-strength, the LN 
pairing-energy also collapses. Therefore, the LN method which is used so extensively in 
the nuclear structure studies appears to be marginally better than the HFB approach and 
deviates substantially from the PHFB results. As is evident from Fig. 1, the PHFB pairing- 
energy is non-zero for all the values of the pairing-strength. In the following, we shall now 
turn to the rotational aspects of the pairing-energy. 

The results of the cranking calculations with monopole-, (monopole plus quadrupole)- 
and delta-interactions are presented in Figs. 2, 3 and 4, respectively. In the upper panel 
of Fig. 2, a marked difference between the PHFB and HFB results are obtained for the 
monopole-interaction at lower rotational frequencies. The PAV on the other hand gives a 
good description for the total-energy and coincides with the PHFB total-energy at lower 
rotational frequencies. However, at higher rotational frequencies it shows a transitional 
behaviour and goes over to the HFB total-energy. The reason for this transition is that 
at higher rotational frequencies, the HFB pairing-energy vanishes and the results of PAV 
become identical to HFB. It should be noted that the mean-field potential, apart from the 
contribution of the pairing- interaction (F-term), is fixed in the present model study and as 
the pairing-energy vanishes the results of projected-energy are same as that of unprojected 
HFB in the PAV method. The LN total-energy is slightly lower than HFB, but is still far 
from the PHFB total-energy. At higher rotational frequencies all the approximate projection 
methods coincide with the PHFB, since the pairing-cncrgy becomes quite weak. 

The HFB pairing-energy, shown in the middle panel of Fig. 2, depicts a transitional 
behaviour as expected. For fiuj < 0.45G, the pairing-energy is finite, but for higher rotational 
frequencies it vanishes. This is a major problem with the HFB theory and is related to 
the fact that it does not consider the particle-number fluctuations which become quite 
important for a mesoscopic system like an atomic nucleus. The PHFB pairing-energy, on 
the other hand, is clearly non-zero for all the rotational frequencies. It does drop at around 
huj — 0.55G and continues to drop with increasing rotational frequency, but decreases very 
smoothly. This smooth drop is due to the fluctuations which are inherent in the projection 
formalism. The PAV pairing-energy is somewhat better at low rotational frequencies, but 
after the HFB phase transition is same as that of HFB. The Lipkin-Nogami pairing-energy 
shows a slightly smoother drop at the transitional frequency, but then goes to zero at a 
slightly higher frequency as compared to HFB. 

The observable quantity which is directly related to the pairing-energy is the aligned 
angular-momentum along the rotational axes (< >)■ The nuclei in the ground-state are 
in a paired phase with maximal pairing-correlations. However, with increasing rotational 
frequency the particles are unpaired and are aligned towards the rotational axes. The 
aligned angular-momentum, therefore, is a measure of the unpaired particles or inversely of 
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pairing-correlations. It is evident from the bottom panel of Fig. 2 that the HFB aligned 
angular-momentum shows a transition from < J.j. >~ to ~ 10. This jump is related to 
the transitional behaviour of the pairing-energy. The results of PAV closely follow the HFB 
results. The LN aligned angular-momentum has a smoother behaviour, but deviates from 
PHFB. The crossing frequency, which is given by the point of inflexion in the slope of the 

< Jx> curve, is same for HFB and PAV. The crossing for LN is slightly higher than HFB, 
but is lower than PHFB by huo = O.IG. 

The results of the projection methods with (monopole plus quadrupolc)-interaction are 
presented in Fig. 3. The overall agreement between the approximate projection methods 
( HFB, PAV and LN ) and the PHFB results is improved as compared to the results with 
a pure monopole-interaction. The results of the total-energy for PAV agree very well with 
PHFB at lower rotational frequencies as in the case of monopole interaction. The LN 
total-energy is slightly improved as compared to HFB, but deviates considerably from the 
PHFB total-energy. The results of the pairing-energy, shown in the middle panel of Fig. 
3, are similar to the results with monopole pairing. The LN pairing-energy is smoother at 
the bandcrossing, but vanishes at higher rotational frequencies. The behaviour of the LN 
aligned angular-momentum, shown in the bottom panel of Fig. 3, is now quite similar to 

< Ja, > of PHFB. In particular, the slope of the curve, which is a measure of the interaction 
between the ground-state band and the s-band, is quite similar to that of HFB, although 
the bandcrossing is slightly lower than that of PHFB. 

The results of the calculations with the full delta-interaction are shown in Fig. 4. The 
total-energy obtained from different approximate particle number-projection methods is con- 
siderably improved as compared to the previous two cases of monopole- and (monopole plus 
quadrupole)-interactions. However, the HFB pairing-energy shows the phase transition from 
finite pairing to zero pairing as in the other two cases. The LN pairing-energy on the other 
hand is quite smooth and is non-zero up to a very high rotational frequency. Due to this 
smooth behaviour, the aligned angular-momentum also shows a smooth behaviour and has 
a similar slope as that of PHFB. 

We have also studied the deformation dependence of the pairing-energy obtained from 
various projection methods. As the deformation increases, the pairing is expected to drop. In 
particular, for superdeformed shapes the pairing is quite weak and the results become quite 
sensitive to the treatment of the pairing-energy. In Fig. 5, the pairing-energy is presented 
as a function of the deformation parameter, n for three different pairing-interactions shown 
in the three panels of the figure. As is clear from the figure, the pairing-energy drops with 
increasing deformation obtained from all the methods. However, the HFB pairing-energy 
goes to zero at around k — 5G with monopole- and (monopole plus quadrupole)-interactions 
shown in the upper two panels of Fig. 5. In the case of full delta-interaction, the HFB 
pairing-energy collapses at a slightly higher deformation value of k, = 5.5G. The LN pairing- 
energy has a smooth behaviour as a function of the deformation, which is quite similar to 
pairing-energy obtained in the PHFB case. However, the LN pairing-energy also goes to zero 
at higher deformation values, wheraes the PHFB pairing-energy is finite for all the studied 
deformation values. 

In Fig. 6, the results of the cranking caculations are presented for a very large deforma- 
tion value of K = 5.5 which is analogous to the deformation of the superdeformed bands. 
The dynamic moment of inertia (J^^-*) shown in Fig. 6 is obtianed from the calculated 
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aligned angular-momentum as a function of the rotational frequency. As is evident from 
Fig. 6 the HFB J^"^^ is quite large as compared to J*^^' obtained from PHFB for all the three 
pairing interactions. The Lipkin-Nogami 

J(2) 

appears to be close to the PHFB J^"^^ at lower 
rotational frequencies. However, at higher rotational frequencies it deviates quite a lot from 
the PHFB J^^^ and as a matter of fact it becomes worse than HFB. 

V. SUMMARY 

In the present study, an attempt has been made to test the quality of the various particle- 
number projection methods by comparing the pairing-correlations and other related quan- 
tities. First of all, it is clear from the present study that the mean-field HFB approach 
appears to be a better approximation with more realistic pairing-interaction. This to our 
knowledge is not discussed in the literature. The HFB results for total-energy and the 
aligned angular-momentum appear to be in better agreement with the PHFB or exact re- 
sults for the delta-interaction. This is in contrast to what is expected from the BCS theory. 
The BCS theory is known to be a better approximation for the simple monopole-interaction. 
This is easy to understand since in BCS theory one only considers the interaction among the 
particles in the time-reversed states, which is appropriate for the monopole-interaction. But 
in a more realistic pairing-interaction, one also has particle-particle correlations which do 
not not necessarily operate in time-reversed states. These extra correlations are obviously 
neglected in the BCS theory. In the HFB theory, on the other hand, all the particle-particle 
correlations are considered in a self-consistent manner. We would also like to point out 
that we have considered all the terms arising from the pairing-interaction. In most of the 
studies, the exchange term arising from the pairing-interaction (F) is not considered since it 
contributes to the already fitted Hartree-Fock potential. This contribution turns out to be 
largest with the delta-interaction. In the absence of this term, the agreement between HFB 
and PHFB is better for the monopole-pairing. 

The usefulness of a projection method is evaluated from its description of the pairing- 
correlations. The mean-field HFB approach suffers from a fundamental problem that it leads 
to a sharp phase transition from a superfiuid to an unpaired phase. This is clear from the 
results presented in the present work, the HFB pairing-energy at critical points goes to zero. 
The major test of a projection method is how well it cures for this transitional problem. 
The PAV method is obviously not meant to solve this transitional problem, since it uses the 
HFB solution to start with. It suffers from the same problems as that of HFB. 

The LN method which is used quite extensively appears to improve the pairing-energy 
somewhat, but is also seen to breakdown in the weak pairing limit. It is observed that as 
a function of the pairing-strength it does have a non-zero pairing-correlations at G = 0.4, 
which is a critical point in the HFB case. However, it should be noted that the pairing- 
correlations derived from LN are quite weak as compared to the PHFB case. As a function of 
the rotational frequency, the LN pairing-energy has a similar behaviour, it does not approach 
zero at the HFB critical frequency, but goes to zero at a slightly higher frequency. It can 
be, therefore, concluded from the present study that Lipkin-Nogami method is marginally 
better than the HFB approach. 

The pairing-energy calculated from the PHFB method is seen to be finite for all the 
studied cases. As already mentioned, the calculated total-energy and the aligned- angular 
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momentum from the PHFB method are in complete agreement with the exact solutions of 
the model Hamiltonian. We consider that the PHFB method is an ideal candidate to replace 
the Lipkin-Nogami method in the nuclear structure studies. It is to be noted that one of the, 
often stated, advantages of the LN method is its computational simplicity |[15|-|T7|. The LN 
method, apart from the evaluation of the A2 Lagrangian parameter, is exactly same as that 
of HFB approach. The PHFB equations have also same structure as that of HFB, except 
that one needs to evaluate the projected-fields. The evaluation of these projected-fields 
requires some extra work. 

We would like to caution that a simple model Hamiltonian has been employed in the 
present work and the question arises whether the conclusions drawn here can be extrapolated 
to a realistic Hamiltonian and model space. First of all, we would like to mention that the 
present model Hamiltonian is not totally unrealistic as it contains most of the essential 
elements of a realistic Hamiltonian. The important difference is in the model configuration 
space. In the present work, we have employed the model space of a single-j shell and 
a realistic model space contains several major shells. Due to this major difference, the 
computation of the projected-fields may take a longer time. We are presently implementing 
the PHFB method in a realistic configuration space and our initial estimate is that the 
projected calculations may increase the numerical work by a factor of three to four as 
compared to the normal HFB case. 
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VI. APPENDIX 



We shall derive here general expressions for particle-number projected averages and their 
variation, valid at zero and finite temperature. The approach is also applicable to other 
projections. Let R be an arbitrary many-body density operator. The average of an operator 
O with respect to R will be denoted as 

{0)n = Tt[OR] , (63) 

where the trace is taken in the full grand canonical ensemble. For a Hamiltonian H satisfying 
[H, N] = 0, the average projected-energy is given by 

E"" = = Jd<p (0) , (64) 

where 

p-i4>N I i4>N\ (Hp'-'l'^\a 
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with J d(f)y^{(j)) = 1. Eq. reduces to (|T8D in the pure limit R = |$)($|. 
Let us consider now a density of the form 



i? = e-^^VTr[e-^^] , (66) 

where K is an hermitian one-body operator. In a fermion system, the most general one-body 
K can be written as 

^ = E + I (/C2°c|c] + /C°fQc,) = ko + \Z^1CZ , (67) 



z-r\ z^-ic^c) jc-r^" ] 

with kf) = |Tr/C^^ and K?^, KP^ antisymmetric matrices. For K hermitian, lO = /C, i.e., 
(IC'Y = (^'')*, = -(/C^O)* = (/C20)t. The norm in (m is 



Tr[e-^^] = e-^^° det[l + e-^^^j^/^ . (68) 

It is convenient to define a generalized density matrix of elements TZmn = {ZlZm)R, that 
comprises both the standard density matrix and the pairing tensors and is a function of /C: 

^-('^ (69) 

where pij = {c]ci)R, Kij = {cjCi)R, kij = {c]cI)r, For K hermitian, TZ^ = TZ, i.e., = p*, 
All averages with respect to R can be expressed in terms of TZ. In particular, 
for an operator O of the form 



{0)r = Tr[0"p + liO^'k + O'^^)] =0o + |Tr[07^] . 
Let us write now the particle-number operator N = Y, 

N = no + lzWZ, Ar=(^J_°J, (70) 

with 2no the dimension of the single particle space. As the product of exponentials of one- 
body operators is the exponential of a one-body-operator, we have e^'^'^ e~^^ = e~^^^'^\ with 
K{(j)) a one-body operator determined by |18[ 



K{<f)) = ko{<p) + lZ^}C{<P)Z, (71) 



/9/C(0) _ i<t>Af-f3JC _ 



g^0A^7^(l _ 7^)-i , (72) 



and /co(0) = ~ icpuo/p. Eqs. (|6^ ) - ([72|) lead then to 

{e^nn = det[l + e-.^]i/^ = det[e"^/^C(0)]V^ , (73) 
C(0) = 1 + (e'<^^- l)7^, (74) 
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which enables the evaluation of Besides, by means of Wick's theorem, can 

be expressed in terms of the elements of the non-hermitian density matrix determined by 

7^(0) ^ ( ^1^] ^ ) = [1 + e^'^^^^r = e^^^7^C-(0) , (75) 

where = {c]ci)R(^^), = {cjCi)^^^), k{(f)) = (c]c|)ij,(<^). In the case (^), 

H{<P) = Tr{[e + |r(0)]p(0) + \[A{<j>)R{<p) + A{<f>)K{<f>)]} , (76) 

with Tij{(f)) = T,k,iVikjiPik{4>), Ajj(0) = ^T,k,iVijkii^ki{<P), Aij{(j)) = lT,k,i'Vkiijiiki{<P)- Let us 
now evaluate the variation of ([6^ ) with respect to TZ. We have 

S{e'^^)R = l{e'^^)RTT{C-\<P)[e''^^ - l]57^} , 
Sn{(f)) = e'^^C~\(j)){6n)C^\(j)), C(0) = l + 7^(e''^•^ - 1) , 
6H{(j)) = |Tr[7^(0)57^(0)] = |Tr[C-l(0)7^(0)e*<^■^C"l(0)57^] , 

where H{(j)) is defined by ^Hmn{(p) = dH{(j))/d7lnm{(p), i-e., 

Eqs. (|77|)-(|78D are valid for any many-body Hamiltonian. In the case (^), h{(j)) = e + r(0) 
and A(0), A(0) coincide with the previous definitions. We finally obtain 

5E^ = lTi[n^Sn] , (79) 

where 7i^, defined by hTi.^^ = dE^ /dTZnm, is given by 



= y d0y^(0)C-i(0)[(e'^^ - l)(i/(0) - E"^) + H(0)e''^^C-^ 



(80) 



For H hermitian, {H^y = . The first term C-^(0)(e^'^^ - l)(i/(0) - ^^) arises from 
the variation of y^{(f)). 

Eq. ([79|) holds for arbitrary variations 57^. For unitary variations 6TZ = i[6yV,TZ], Eq. 
([79|) leads to 6E'^ = — |iTr{[7Y^, 7?.] SW}. Stabihty against these variations leads then to 
the necessary condition 

[7^^,7^]=0, (81) 

which has the same form as the standard HFB condition [Q] except that the HFB matrix Ti, 
has been replaced by Ti.^ . 
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All previous expressions hold actually in any basis of single particle or quasiparticle 
states, after appropriate transformation of M . The matrix /C can be diagonalized with a 
transformation (|), i.e., Z' = W^Z, with {Z')^ = (a^a). We may then write JC = W/C'>V^ 
with )C' diagonal, and 

n = wR'w\ ^' = 1 - / ) ' ^^^^ 

with /fc = {alak)R = (1 + e^^'=)~^ and E^. the positive eigenvalues of /C. 

T = Urmt. For T = P'^ ^ 0, fk ^ (assuming Ek ^ 0) and R |0)(0|, with |0) the 
vacuum of the quasiparticles a^. In this limit, 

^' ^ ( 1 ) ' ^ ^ ( t/yt uw) ' (^^^ 

implying p = (VV^)*, k = {VW)* = k\ and 7^^ = TZ, i.e., + kk = p, pn = Kp^ . Using 
these relations, we obtain the simplified expressions 

det[e*^/2(3(0)]^/2 ^ det[e*^D(0)]^/2 ^ ^ ^i^^ ^ ^-i4>^i _ ^) ^ 

l<-[<P) ^ l^e-i<^^/^-i(0) i-e^^[D-\<P)pY ) ' ^ ^ 

withC-^(0) = l+(e-^<^•^-l)7^(0), C-^(0) = l-e-^<^•^7^(0)(e*'^•^-l). With the replacements 
C(0) = e^<^D-i(0), k(0) = Kt(0)^ A(0) = At(0), Eq. (||) leads to Eqs. (|6D-(||). 

For T = variations that preserve the idempotent condition TZ^ = TZ, we may employ a 
simplified matrix 7i^, which can be directly derived from the above expressions. Its elements 
are 

= J d(j) {(f)) D-\(f)){i sin (j){H{(j)) - E^) 

+ [h{(f)) - i sin (f){e''''KA{(j)) + e-'^A{(j))R)]D-\(j))} , (85) 
- ^rf0y^(0)e-*^D-i(0)A(0) , = / rf0y^(0)e*^A(0)D-i(0) . 



In the case of Hamiltonian (|T|), with the previous replacements these expressions lead to 
Eqs. dl) ^ 
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FIGURES 



FIG. 1. The results of the total energy (Efot) and the pair-energy {Epair) as a function of the 
strength of the pairing-interaction (G) using various particle-number projection methods. 

FIG. 2. The results of the total energy (Etot), the pair-energy (Epair) and the alignment (Jx) 
for six-particlcs in a deformed j = 11/2 orbitial using the monopole interaction. The results are 
obtained with Hartree-Fock-Bogoliubov (HFB), projection after variation (PAV), Lipkin-Nogami 
(LN) and projected HFB (PHFB) approaches. 

FIG. 3. The results of the total energy (Etot), the pair-energy (Epair) and the alignment (Jx) 
using the monopole plus quadrupole pairing interaction. 

FIG. 4. The results of the total energy (Etot), the pair-energy (Epair) and the alignment (Jx) 
using the full delta- interaction. 

FIG. 5. The results of the function of the deformation parameter 

(k). (a) gives the results with monopole-pairing, (b) presents the results with (monopole plus 
quadrupole)-interaction and (c) gives the results with full delta-interaction. 

FIG. 6. The dynamic moment of ineratia for k = 5.5. (a) gives the results with 
monopole-pairing, (b) presents the results with (monopole plus quadrupole)-interaction and (c) 
gives the results with full delta-interaction. 
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